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Quantum Valence Criticality as Origin of Unconventional Critical Phenomena 
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It is shown that unconventional critical phenomena commonly observed in paramagnetic metals 
YbRli2Si2, YbRli2(Sio.95Geo.o5)2, and /3-YbAlB4 is naturally explained by the quantum criticality of 
Yb- valence fluctuations. We construct the mode coupling theory taking account of local correlation 
effects of f electrons and find that unconventional criticality is caused by the locality of the valence 
fluctuation mode. We show that measured low-temperature anomalies such as divergence of uniform 
spin susceptibility \ ~ with £ ~ 0.6 giving rise to a huge enhancement of the Wilson ratio and 
the emergence of T-linear resistivity are explained in a unified way. 

PACS numbers: 71.27. +a, 75.20.Hr, 71.10.-w, 71.20.Eh 



Quantum critical phenomena have been one of the cen- 
tral issues in condensed matter physics for the past two 
decades. The nature of the quantum critical point (QCP) 
emerging when the magnetically-ordered temperature is 
suppressed to absolute zero has been intensively studied 
and the role of spin fluctuations in the critical phenomena 
has been well understood 0-0]. 

However, anomalous critical phenomena, which do not 
follow the conventional spin-fluctuation theories 0-0 
have been discovered in paramagnetic metal phases in 
YbRh 2 Si 2 0, 0, YbRh 2 (Si . 95 Geo.o 5 )2 0, |, and (3- 
YbAlB-4 @. The most striking anomaly is that low- 
temperature uniform spin susceptibility exhibits diver- 
gent behavior x{T) ~ with the anomalous criti- 
cal exponent £ = 0.6 in YbRh.2(Sio. 95660.05)2 and £ = 
0.5 in /3-YbAlB4 in spite of no sign of a ferromag- 
netic phase nearby. In these materials, the Sommer- 
feld constant exhibits the logarithmic divergence j c = 
C c /T ~ -InT 0, i, |H, giving rise to a large Wil- 
son ratio, i.e., a dimensionless ratio of x to j e , Rw = 
17.5 in YbRh 2 (Sio.95Ge .o5)2 0, and R w = 6.5 in [3- 
YbAlB>4 0], exceeding the conventional strong-coupling 
value i?w = 2. The linear-T dependence of low-T re- 
sistivity emerges in a wide-T range in YbRh 2 Si 2 0, 
YbRh 2 (Si .95Ge .o5)2 0, and also /3-YbAlB 4 0]. These 
observations suggest that there exists a new class of ma- 
terials showing similar critical phenomena, which are un- 
conventional. 

So far, to explain YbRh 2 Si 2 and YbRh 2 (Sio. 9 5Geo.o5)2, 
theoretical efforts have been made [Iol - [l2j |. In particular, 
a scenario asserting that f electrons undergo a localized 
to itinerant transition was extensively discussed [l(| E3 ■ 
However, YbRh 2 Si 2 shows a large Sommerfeld constant 
as ~ 1.7/JmolK 2 even inside of the antiferromagnetic 
(AF) phase 8]. This fact indicates that heavy quasi- 
particles are responsible for the AF state. Indeed, a 
band-structure calculation showed evidence contrary to 
the scenario by demonstrating that a tiny valence change 
of Yb can exp lain the Hall-coefficient measurement in 
YbRh 2 Si 2 jL3|, which seems to be a basis of the above 
scenario jl4j . 

Recently, in /3-YbAlB-4 the valence of Yb has been de- 



tected as Yb +2 75 (0.75 4f hole per Yb]_at T 
suggesting strong valence fluctuations 
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20 K, 
Close rela- 
tion of anomalous critical phenomena and valence fluc- 
tuations has been also indicated in Ceo.g-ajTho.iLaa; [l6[ 
and in YbAuCu4 [TtJ ■ Both are related to typical valence- 
transition materials: One is Ce metal, well known as y-a 
transition [l8| , and the other is YbInCu4 [l9[ , both show- 
ing a discontinuous valence jump of a Ce and Yb ion, 
respectively, when T and P are varied. Since the first- 
order valence transition is an isostructural transition, the 
critical end point exists in the T-P (and -chemical com- 
position) phase diagram, as in the liquid-gas transition. 

At x ~ 0.1 in Ceo.g-zTho.iLaa,, at which the critical 
end point is most suppressed and is close to T = K [l(| , 
critical phenomena rising from the quantum critical end 
point were revealed: T-linear resistivity emerges promi- 
nently and uniform spin susceptibility is enhanced at 
low T giving rise to a large Wilson ratio, i?w ~ 3. In 
YbAuCu4 , the uniform spin susceptibility is enhanced as 
x(T) ~ 7 1 — 0-6 as t decreases in spite of the fact that 
the AF transition takes place at Tn = 0.8 K [20j, sim- 
ilarly to YbRh 2 Si 2 . Furthermore, the sharp Yb- valence 
crossover temperature T*{H) is induced by applying a 
magnetic field [17(, suggesting that YbAuCu4 is located 
in the vicinity of the quantum critical end point of the 
valence transition. It should also be noted that the T-H 
phase diagram of YbAuCu4 pjj closely resembles that 
of YbRh 2 Si 2 0. The H dependence of the crossover 
temperature T*(H), emerging in several physical quan- 
tities 0, whose origin is unclear in YbRh 2 Si 2 is quite 
similar to the T*(H) in YbAuCu 4 . 

These observations strongly suggest the importance 
of quantum criticality of valence transition as a key 
mechanism of unconventional critical phenomena. From 
this viewpoint, in this Letter, we resolve this outstand- 
ing puzzle by showing that (a) uniform spin suscepti- 
bility diverges with anomalous criticality x{T) ~ 
with 0.5 < C < 0.7 in paramagnetic metals even with- 
out proximity to a ferromagnetic phase and (b) T-linear 
resistivity emerges in the wide-T range. 

Let us start our discussion by introducing a minimal 
model which describes the essential part of the Ce- and 
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Yb-based systems in the standard notation: 



H = H r + H f + H ] 



hyb 



H, 



(1) 



where H c = Ek<x £ k4 CT c k CT , H i = £ f E«x ™Lx + 

and i?(7 fc = t/f c E*Li "i n i ■ The C/f c term is the Coulomb 
repulsion between 4f and conduction electrons (holes) in 
Ce (Yb) systems, which is considered to play an impor- 
tant role in the valence transition [2l| . 

To consider correlation effects by Ug, we emp loy the 
slave-boson large- N expansion framework J22|. The 
Hamiltonian Eq. (JXJ) is generalized to the case with TV-fold 
degeneracy from a = J [, 4,, and the slave-boson operator hi 
is introduced to eliminate the doubly-occupied state for 
Us — > 00 under the constraint n\ T J Y '•' 
Lagrangian is written as C = Co + C 



„„ , A/,,7,, = 1. The 
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where Ak is the Lagrange multiplier to impose the con- 
straint, and £k = £k + an< ^ ^L-k' = ( £ f + if 1 ) ^kk' + 

^T^Ak-k'- We here separate £ as Cq and £' to perform 
the expansion with respect to the U{ c term after taking 
account of the local correlation of the Us term. 

For exp(— Sq) with the action Sq — cItCq(t), the 
saddle point solution is obtained via the stationary con- 
dition SSq — by approximating spatially uniform and 
time independent ones, i.e., A q = A<5 q and 6 q = 6<5 q . 
The solution is obtained by solving mean-field equations 
dSo /dX — and dSo jdb — self-consistently. 

For S' = J^drCir), we introduce the iden- 
tity applied by a Stratonovich-Hubbard transfor- 
mation e _s ' = fDfexp[J2 im Jo dr{-^tp im (T) 2 + 
i^(cimfi m ~ fim4m)} ( Pim( T ))]- Tne partition function 
is expressed as Z = J V(cc' i //V) exp(— 5) with S = So+ 
S'. By performing Grassmann number integrals for cc^ 
and fp, we obtain Z = J T><pexp(—S[<p]) with S[ip] — 
E m E ff Vm(«)^(-g) - TrlnhGo- 1 + V] - ^A|6| 2 + 
/S/VqoV^sA. Here, the abbreviation q = (q, iwj) 
lui = 21ttT is used, and Go and V are defined as 



with 
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with Vmkk' = V pNN s V m ( k ~ k ')' respectively. Here, k = 
(k, is n ) with £„ = (2n + 1)ttT, and each matrix element 
of Go is defined as (G )n = Gfj c , (G ) 12 = Gfj f , (G ) 2 i = 



G C , and (G ) 22 = GJf. By using Trln[— Gq + V] 
Trlnl-Gg 1 ] - £~ =1 ^Tr(G f we obtain 
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i(qi)<fm(q2)<Pm(q3)Vm(q<l)$ E* + 



(2) 



Here, constant terms independent of (p m {q) are omitted 
in Eq. ([2|) since they merely shift the origin of the free 
energy of the system. The coefficient of the quadratic 
term is given by 



2E/f c 
N 



{xg^twO-X^fatwi)} 



,(3) 



where (^ iuj i) = ~ w s E k ,n G o P ( k + Q, + 

iwj)GQ 5 (k, ie n ). Since long wave length |q| -C 9 C around 
q = and low frequency |u>| lo c regions play dominant 
roles in critical phenomena with q c and w c being cutoffs 
for momentum and frequency in the order of inverse of 
the lattice constant and the effective Fermi energy, re- 
spectively, Qi for £ = 2,3, and 4 are expanded for q and 
cj around (0, 0): 



yOtj3 I 



ft 2 (q,iw ; ) ^r) + Aq 2 + C 
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(4) 



where rj = [7 fc [l - ^{^-(0,0) - xf f (O,0)}], 

^3(91,92,93) ~ Va/^/3N B , and f2 4 (9i, 92, 93, 94) ~ 
va/(PN„). 

Different from ordinary critical phenomena of spin fluc- 
tuations [H-Q, there appears a cubic term in Eq. in 
general for the valence fluctuation case (23|. Let us here 
apply the Hertz's renormalization-group procedure Q to 
S[tp]: (a) Integrating out high momentum and frequency 
parts for q c /s < q < q c and uj c /s z < 10 < u) c , respec- 
tively, with s being a dimensionless scaling parameter 
(s > 1) and z the dynamical exponent, (b) Scaling 
of q and w by q' — sq and u)' = s z u>. (c) Re-scaling 
of ip by (p'(q' = s a ip(q l /s,u)' /s). Then, we derter- 
mined that to make the Gaussian term in Eq. © scale 
invariant, a must satisfy a = —(d + z + 2)/2 with d spa- 
tial dimension and the dynamical exponent z = 3. The 
renormalization-group equations for coupling constants 
Vj are derived as = [6 — (d + z)] v% + 0(v%), and 
^ = [4 — (d + z)} Vi + 0(v\), for cubic and quadratic 
terms, respectively. By solving these equations, it is 
shown that higher order terms than the Gaussian term 
are irrelevant 



lim Vj(s) = for j > 3 



(5) 
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for d + z > 6. For the case of d = 3 and z = 3, 
it is shown that the cubic term is marginally irrele- 
vant [2l|. Hence, the universality class of the critical- 
ity of valence fluctuations belongs to the Gaussian fixed 
point. This implies that critical valence fluctuations are 
qualitatively described by the RPA framework with re- 
spect to £/f c . The coefficient of the Gaussian term in 
Eq. @ is nothing but the inverse of the valence suscep- 
tibility f2 2 (q, i^/) = Xv(q, iw/,) -1 . Since evaluation of 
Xo° c (q, ioJi) and Xo fcf (q ; * w i) using the saddle point solu- 
tion for exp(— S ) concludes Xo cc ^ Xo fcf ( see Fig. 2 and 
text below) , it turns out that Xv is expressed by the RPA 
form Xv(q,*wi) = dT(T T n f (q,T)n f (~q,0))e l ^ T « 
Ufc'll - ^x£ cc (q,^)r\ as shown in Fig.fU 
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FIG. 2: (color) (a) q dependence of Xo cc (9>0) (solid line) and 
Xo fcf (9, 0) (dashed line) calculated by using saddle point so- 
lution of exp(— So) for gf = —1 (green), —0.5 (red), and 0.0 
(black), (b) fit vs £f. (a) and (b) are results for D = 1, 
V = 0.5, and t/ ff = oo at n = 7/8. 
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FIG. 1: Feynman diagrams for dynamical valence suscepti- 
bility and dynamical spin susceptibility for f electrons. Solid 
lines and dashed lines represent the f- and conduction-electron 
Green functions, Go and G§ c , respectively. Wiggly lines rep- 
resent [7f c . 

An important consequence of this result is 
that dynamical f-spin susceptibility Xr (q> ) = 
f£ dr(r T S* f + (q,T)5 f ~(-q,0))e i " ;T has a common struc- 
ture with Xv i n the RPA framework as shown in 
Fig. [TJ At the quantum critical end point of the valence 
transition, namely, the QCP, the valence susceptibility 
Xv(0, 0) diverges. The common structure indicates that 
Xf (0,0) also diverges at the QCP. The uniform spin 
susceptibility is given by x ~ xl ~ l^iSf Xf (0) 0) with 
xl uniform f-spin susceptibility, /xb the Bohr magneton, 
and c/f Lande's g factor for f electrons. This gives a 
qualitative explanation for the fact that uniform spin 
susceptibility diverges at the QCP of valence transi- 
tion under a magnetic field, which was shown by the 
slave-boson mean-field theory applied to Eq. (UJ) [24j . 
Numerical calculations for Eq. {TJ in d = 1 by the 
DMRG [24| and in d — 00 by the DMFT [if also 
showed the simultaneous divergence of Xv and uniform 
spin susceptibility under the magnetic field, reinforcing 
the above argument based on RPA. 

The other important point of the present theory is that 
the "unperturbed" term Lq, i.e., Go, already contains 
the local correlation effect by Us- This effect plays a key 
role in critical phenomena in Ce- and Yb-based systems, 
which will be shown below to be the origin of the uncon- 
ventional criticality. The local correlation effect emerges 
as dispersionless, almost g-independent Xo° c ( ( l'0) an d 
Xo fcf (q, 0) in Eq. (0), as shown in Fig. [2fa). Here, the 
saddle point solution for exp(— So) is employed for a typ- 
ical parameter set of heavy-electron systems: D = 1, 
V = 0.5, and Us = 00 at total filling n = 7/8 with 
£k = k 2 / (2mo) — D and n = fif + n c , where fif and n c are 



the number of f electrons and conduction electrons per 
"spin" and site, respectively. The bare mass mo is chosen 
such that the integration from — D to D of the density 
of states of conduction electrons per "spin" is equal to 1. 
This local nature is reflected in the inverse of valence sus- 
ceptibility in Eq. ((4]) as an extremely small coefficient A. 
We note here that this flat-g result is obtained not only 
for deep £f with fif = 1/2 in the Kondo regime, but also 
for shallow Ef with fif < 1/2 in the valence-fluctuating 
regime (see Fig. [2Jb)). Here, we note that the c-f hy- 
bridization is always finite. 

To clarify how this local nature causes unconventional 
criticality, we construct a self-consistent renormaliza- 
tion (SCR) theory for valance fluctuations. Although 
higher order terms Vj (j > 3) in S[tp] are irrelevant 
as shown in Eq. ([5]), the effect of their mode couplings 
affects low-T physical quantities significantly as is well 
known in spin-fluctuation theories [lHj]. To construct 
the action using the best Gaussian taking account of the 
mode-coupling effects up to the 4-th order (j < 4) in 
S[tp], we employ Feynman's inequality on the free en- 
ergy: F < F cS + T(S - S cS ) cS = F( V ). Here, S e ffM = 
\ E m E q ,i( ? 7 + yll ? 2 + G g |a; ; |)|^ m (q, kj/)| 2 , and 77 is deter- 
mined to make F(j]) be optimum. By optimal condition 
dF(r])/dr] = 0, the self-consistent equation for 77, i.e., 
the SCR equation, is obtained: 77 = 770 + 3«4(iy9 m ) 2 ff /N s , 
where ((p 2 m ) eS = rE q j(l + Aq 2 + C q \ui\)- 1 . Here, we 
write (<p m ) e ff in a general form using C q , which is given by 
C q = G j max{g, l^ 1 } with l[ being the mean free path by 
impurity scattering [26]. When the system is clean, i.e., 
C q — C/q, the SCR equation in d — 3 in the Ag B < 7/ 
regime with q& being the momentum at the Brillouin 
Zone is given by 



y = yo + ^y\t 



62/ 2?7 Jo 



jt 

6y _ 



(6) 



where y = v/(Mb), t = T/T , T = Aql/pnC), x = 
q/qB, x c = q c /qB, and 7/0 and yi are constants. When 
y^>t,y<x t 2 / 3 is obtained from Eq. ^ at the QCP with 
7/0 = 0. This indicates that the valence susceptibility 
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FIG. 3: (color online) (a) Numerical solutions of Eq. ([6]) for 
yo = 0.0 (at QCP), 0.3, and 0.6 at yi = 1 and x c = 1. 
(b) Electrical resistivity p(T) calculated by using y(t) in (a). 
Dashed line represents the linear-t fit. 

shows unconventional criticality Xv(0, 0) = rf 1 oc t~ 2 ^ 3 . 
Figure shows numerical solutions of Eq. ^ . Note 
here that the coefficient A is quite small as shown above, 
giving rise to quite small T (= Aq^/(2nC) <C Tp) so 
that the region of t = T/T$ shown in Fig. [3] corresponds 
to that of T <C Tp ~ O(D). Hence, a wide range of 
t = T/To is shown in the plot. The least square fit of 
the data for 5 < t < 100 gives y oc f - 551 . Since the 
Gaussian fixed point ensures the simultaneous divergence 
of valence susceptibility and uniform f-spin susceptibility 
as discussed above, xl shows divergent behavior ~ 
with 0.5 < C < 0.7 depending on the temperature range 

in agreement with experiments in YbRh2(Sio.95Geo.os)2 
and /3-YbAlB 4 . We note that the NMR or NQR re- 
laxation rate is shown to be (TiT)^ 1 ~ Xs(t) °c t~^, 
which also quantitatively agrees with (TiT) -1 ~ y-OS 
in YbRh 2 Si 2 H3]. 

We should note here that in the T — > limit, although 
it may be experimentally difficult to access such a low 
temperature overcoming the smallness of A, the SCR 
equation follows the conventional z = 3 type [l|, H, ||J 
in the clean system, giving rise to y oc i 4 / 3 . Then, at the 
QCP (yo = 0), as t decreases, a crossover from y oc i 2 / 3 
to y oc i 4 / 3 occurs. In reality, however, because of the 
smallness of A, the low-T range is extremely elongated 
by the relation t = T/Tq with Tq oc A, which makes it 
possible that unconventional criticality dominates over 



the experimentally accessible low temperature region. 

We note that the electrical resistivity p(T) shows a 
T- linear dependence in the regime t > 5 (y > 1) where 

Eq. §§§ is applicable, as shown in Fig. [3fb). Here, 
following a formalism of Ref. (28[, p(T) is calculated 
as p(T) oc i dujuin(uj)[n(uj) + 1] J^ c dqq 3 Imx^(q,uj) 
withn(o;) = l/(e' 3LJ — 1) being the Bose distribution func- 
tion, and x^iQiW) = (r) + Aq 2 — iCqUj)^ 1 , a retarded va- 
lence susceptibility. Here, y(t) in Fig.^a) is used for the 
clean system C q — C/q, and the normalization constant 
is taken as 1 in the p(t) plot. The emergence of pit) oc t 
behavior can be understood from the locality of valence 
fluctuations: In the system with a small coefficient A, 
where the local character is strong, the dynamical expo- 
nent may be regarded as z = oo when we write C q in 
a general form as C q = C/q z ~ 2 . By using this expres- 
sion in Xvfe '-') m the calculation of piT) for z = oo, 
we obtain p(T) oc T toward T — > K. This result indi- 
cates that the locality of valence fluctuations causes the 
T-linear resistivity. The emergence of p(T) oc T by va- 
lence fluctuations was shown theoretically on the basis of 
the valence susceptibility \v which has an approximated 
form for z = oo in Ref. 29] . Our present formulation and 
the renormalization analysis provide a reasonable grou nd 
for the Xv introduced phenomenologically in Ref. [29(. 

The evaluation of the quasiparticle self energy for a va- 
lence fluctuation exchange process by using the Xv shows 
that ReS(e) oc eln(e), which leads to a logarithmic- 
T dependence in the specific heat C/T for a certain-T 
range [30|. The detailed T dependence of C/T will be 
discussed in a separated paper. 

In summary, we have shown that unconven- 
tional criticality commonly observed in YbRh-2Si2, 
YbRh 2 (Sio.9 5 Geo.o5)2, /3-YbAlB 4 , (V. ;1! , , Th,,.; I.;,,.. and 
YbAuCu4 is naturally explained from the viewpoint of 
the quantum valence criticality. The locality of the Yb 
or Ce valence fluctuation mode rising from local correla- 
tions of the 4f electrons is revealed to be the key origin. 
It is noted that our results may be regarded as an explicit 
manifestation of the marginal Fermi liquid theory [3l| . 
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